Manipulation of qubits in non-orthogonal collective storage modes 
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We present an analysis of transfer of quantum information between the collective spin degrees of 
freedom of a large ensemble of two- level systems and a single central qubit. The coupling between the 
central qubit and the individual ensemble members may be varied and thus provides access to more 
than a single storage mode. Means to store and manipulate several independent qubits are derived 
for the case where the variation in coupling strengths does not allow addressing of orthogonal modes 
of the ensemble. While our procedures and analysis may apply to a number of different physical 
systems, for concreteness, we study the transfer of quantum states between a single electron spin 
and an ensemble of nuclear spins in a quantum dot. 
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I. INTRODUCTION 

Storage of quantum states in the collective quantum 
degrees of freedom of a large ensemble of identical two- 
level systems combines the advantage of the long coher- 
ence lifetime of microscopic systems with the strong cou- 
pling to auxiliary quantum degrees of freedom due to col- 
lective enhancement. Numerous implementations have 
been studied, from storage of optical states of light in 
optically dense atomic ensembles [1-3] to the coupling of 
single electronic spin states with nuclear spin ensembles 
[4-7], and transfer of superconducting qubit states via 
quantized cavity fields to rotational states in molecular 
ensembles [8], and to collective states in electronic and 
nuclear spin states [9-12]. 

For applications in quantum computing and communi- 
cation one needs the capability to store and manipulate 
several qubits. While this can be achieved by applica- 
tion of a separate ensemble for each qubit, there has also 
been a number of proposals to identify orthogonal col- 
lective excitation modes in a single ensemble systems for 
multimode storage and manipulation. 

Multimode storage of light has thus been demonstrated 
in inhomogeneously broadened media [13, 14], while use 
of different molecular [15] and atomic [16, 17] states can 
be used to collectively store separate qubit or oscillator 
states. With spatially extended media, one has the pos- 
sibility to apply ideas from holography and store exci- 
tation patterns with different spatial periodicities [9, 18] 
and also to use the dephasing caused by inhomogeneities 
to address different collective spin superposition states 
[19, 20]. The weak coupling of nuclear spins to their 
surrounding host material and to external perturbations 
make them ideal candidates for long time storage, but 
the same weak coupling makes it difficult to establish 
and address an independent set of collective nuclear spin 



modes, unless, as in [20], one may use the phase evolution 
of an electron spin ensemble to develop collective super- 
position phases and only subsequently transfer them to 
the nuclear spin ensemble. 

In this paper, we assume that some controllable inho- 
mogeneity is available in the coupling strengths of our 
central qubit to the ensemble members, but this inho- 
mogeneity is insufficient to switch the coupling between 
orthogonal and independent collective excitation degrees 
of freedom. We present a method that allows us to ef- 
fectively define and address independent qubits in such 
a system. The analysis may apply for a variety of phys- 
ical systems, but for concreteness we consider a single 
electron, captured in a generic semiconductor quantum 
dot, and coupled through the hyperfme magnetic dipole 
interaction to the nuclear spins in the bulk of the quan- 
tum dot [4]. We assume that the spatial wavefunction of 
the electron can be manipulated by external fields, and 
since the spin-spin coupling depends on the electron den- 
sity at the site of every nucleus, we can hence manipulate 
the precise form of the collective coupling to the nuclear 
spins. 

In Sec. II, we present the physical model and we intro- 
duce the notation and concepts used in the manuscript. 
In Sec. Ill, we identify different, but non-orthogonal nu- 
clear spin modes, that couple to the electron spin. In Sec. 
IV, we present a procedure that allows addressing of two 
orthogonal modes by suitably timed sequences of interac- 
tion with the non-orthogonal spin modes. In Sec. V, we 
investigate the fidelity of our protocol, and we propose 
a high probability heralding scheme which significantly 
improves the fidelity. In Sec. VI we conclude and discuss 
our results. 



II. AN ELECTRON SPIN IN A NUCLEAR SPIN 
BATH 
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We consider the situation depicted in FIG. 1, of the 
spin degree of freedom associated with a spatially con- 
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fined electron wavefunction. The spin of the electron in- 
teracts with external magnetic fields and with the mag- 
netic field created by the nuclear spin ensemble within 
the spatial extent of the electronic state. 




FIG. 1. Quantum dot realized in the tight potential well in 
the plane of a hetero junction of two semiconductor materials 
with different bandgaps. In a lightly doped or even intrinsic 
material with very few impurities, the mobility of the trapped 
electrons can be high and transverse confinement within the 
plane can be established by the potential from electrodes on 
the surface of the semiconductor. The figure indicates that 
the electron wave function, and hence, the interaction between 
the electron and nuclear spins are spatially dependent. 

The magnetic dipole-interaction between the electron 
spin S and the ith nuclear spin Ij is given by 

Hi = giS ■ li 

= 9i{l(S + lr + S-lt)+S,I,,^, (1) 

where the interaction strength, gi, depends on the geom- 
etry of the quantum dot. In particular cx |^?n(rj)| , 
where tp g (ri) is the spatial wave function of the electron 
at the position of the ith nucleus. We assume, that the N 
nuclear spins are perfectly polarized in the — z-direction, 
i.e. to/ i = —Iq. In that limit the total Hamiltonian, 
H g = J^. Hi , can be written as 



Hg = V -^r~ 9 ( s+h ° + s ~ b t) + • b °h- (2) 

In Eq. (2) we have implemented the Holstein-Primakoff- 
approximation, which describes the collective nuclear 
spin excitation by bosonic creation and annihilation op- 
erators. This approximation assumes that the number of 
excitations is substantially less than the total number of 
spins, and as we will consider only 0,1 and 2 excitations 
within an ensemble of thousands of spins, it is perfectly 
valid. In Eq. (2), g = (1/N) ^ \gi\ 2 , is the rms- value 
of the coupling strengths {gi}, 7 is the effective Lande- 
factor of the electron, which depends on the quantum dot 
geometry and host material, and the Overhauser-field, 
Boh, is given by the expression 

1 N 

Boh = ^ 9iIz,iZ- (3) 

With the above definitions we can write the creation 



and annihilation operators explicitly as 

1 N * 1 N 

bt = ^=Y,^<rt and 6fl = J= V) (4) 

Let us now consider the effect of these operators in de- 
tail. We have assumed perfect polarization of the nu- 
clear spins, so we can define our ground state as |0) = 

|0 1 2 ...0 JV ) = |(-Jo)i(-io)2...(-/o)jv>, while, for a 
state with a single excited spin, we write, | . . . lj . . .) = 

|...(-/o + l)<.-.). 
If we apply the Holstein-Primakoff creation operator 

to the ground state we get 
1 N 

|l fl ) = 6t| )= £^| 0l 2 ...l i ...0 JV ), (5) 

V a i=1 g 

which can be interpreted as a single collective excitation 
in the nuclear spin-ensemble. This excited state is char- 
acterized by the set {gi} that describes the pattern of the 
amplitudes, with which the individual nuclear spins have 
been flipped. 

If we inspect the Hamiltonian in Eq. (2), we see that 
it consists of two different terms. The first term is a 
flip-flop-interaction that conserves the total spin projec- 
tion. With the Holstein-Primakoff approximation for 
the nuclear spin states, this interaction is equivalent 
to the Jaynes-Cummings Hamiltonian, effectively cou- 
pling a two-level system (the electron spin) to a har- 
monic oscillator (the excitations in the nuclear spin- 
ensemble). The second term in the Hamiltonian de- 
scribes the electron spin precessing in a magnetic field, 
and we will include also an external field, such that 
Boh -> Boh + B oxt = B off . By adjusting the exter- 
nal field the flip-flop-interaction can be turned on and off 
resonance. 



III. MULTIPLE OSCILLATORS 

As proposed in [4], the nuclear spin degree of freedom 
has a very long lifetime and therefore the quantum in- 
formation represented by the electron spin may be trans- 
ferred to the nuclear ensemble for robust long-time stor- 
age. In this manuscript our goal is to store more than 
a single qubit in the same nuclear spin ensemble. This 
is possible if we can address orthogonal collective spin 
wave modes, as e.g., done in [9, 18, 19]. Here, however, 
the effectively coupled nuclear spin ensemble is confined 
to the spatial volume occupied by a single electron, and 
we do not have the same means to address orthogonal 
plane wave modes. 

We may, however, perturb the spatial wavefunction of 
the electron by applying an electric field with a compo- 
nent parallel to the plane of the quantum dot or we may 
excite the electron to another motional bound state in 
the quantum dot, see also FIG. 2. In this altered state, 
the wavefunction is modified, 4>g{ T ) ~^ (r), and since 
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the coupling strengths between the electron spin and the 
individual nuclear spins depend on the electron spatial 
probability distribution the electron spin hence couples 
to a different collective spin degree of freedom character- 
ized by the coupling strengths, {hi} oc \iph{*i)\ 2 - The 
Hamiltonian which governs this new interaction, Hh, is 
otherwise analogous to Eq. (2), but with the creation and 
annihilation operators replaced by the ones correspond- 
ing to the new oscillator. 




x 



FIG. 2. This figure illustrates how the electronic wavefunction 
of an electron in the ground state (dashed line) of a quantum 
dot can be modified (solid line) by the application of an elec- 
tric field (a), and by excitation of the electron (b). 

By controlling the electron spatial wavefunction we 
thus have a choice between two different storage modes, 
but unless the modes can be addressed independently of 
each other, this does not in itself provide the ability to 
use the nuclear spin ensemble cts ct rc gister for the storage 
of two different qubits. The requirement of independence 
is met if the mode operators for the two spin modes com- 
mute, and in particular, 

'b g ,b[]=0. (6) 

Using the Holstein-Primakoff-approximation, we obtain 

r i 1 N 

which in turn is equal to the overlap (l g |l/,) between the 
single excitation quantum states, |1 Q ) = &J,|0), ot = g, h. 

The states |l g ), |1^) span a two dimensional Hilbert 
space, and it is possible to construct a superposition of 
the states and that is orthogonal to see 
FIG. 3. We label this state and using the Gram- 
Schmidt-method we find an explicit form of \l g )'- 

\l g ) = (l-{l g \l h ) 2 )- 1 / 2 (\l h )-{l g \l h )\l g )^ (g) 

The overlap = ensures that the operator b- 

creating from the fully polarized spin state, indeed 
commutes with the collective raising and lowering oper- 
ators (4) for our original spin wave mode, 



and it is possible to exchange quantum states between 
the electron spin and the g-oscillator without modifying 
the state of the ^-oscillator. 




FIG. 3. A geometric illustration of the mode orthogonal- 
ity problem. We can selectively address one of two different 
oscillators, g- and h. Their first excited states span a two- 
dimensional Hilbert-space but they are not orthogonal. We 
identify a scheme to construct and address a linear superpo- 
sition of the two modes, g, that is orthogonal to g. 



IV. ADRESSING THE QUANTUM MEMORY 

We will use the (/-oscillator mode as the read-in/read- 
out-mode for direct transfer of qubit states between the 
electron spin and the nuclear spin ensemble. When a 
quantum state has been transferred to the g-oscillator, 
subsequent transfer of the state into the g-oscillator can 
be accomplished with the repeated application of the 
pulse sequence 

using precisely the Hamiltonian interaction operators 
available to us, when the electron occupies its two possi- 
ble spatial states. 

The switching between the operators is obtained by 
changing the electronic motional state between ip g and 
iph- By adjusting the external magnetic field the flip- 
flop-interaction can be tuned off resonance while the elec- 
tronic state is manipulated, and when the interaction is 
off resonance the electron spin will precess around the 
z-axis. If we let the electron precess tt radians the S x - 
and S^-operators, and therefore also the S + - and S~- 
operators, will reverse sign. When we return on reso- 
nance, we can thus also obtain the change of sign of the 
interaction Hamiltonian (2) as needed in every second 
application of H g and Hh in (10). 

If we let t be a small time interval, we can use the 
Campbell-Baker-Hausdorff-formula to rewrite the opera- 
tor and, ignoring terms of order higher than r 2 , we get 

Ur ~ e -W*,B,]r' t (u) 

If we apply the sequence N times, the resulting time 
evolution operator is 

U N = {U T ) N = e- N ^ H ^ 2 = e- mT (12) 

where 

R= [H ]}1 H g \ L ^ T ^ mT ^ (13) 
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This means that the evolution corresponds to that caused 
by the effective hamiltonian H in a time interval T, which 
is just defined as the total time of the sequence. 

With the resonant flip-flop-interaction the two Hamil- 
tonians write 

H g = n g (s+b g + s-6t) , H h = n h (s+b h + s-b{) , 

(14) 

and we can evaluate the commutator in Eq. (f3) and 
calculate H. Defining the annihilation operator b g — 
(1 - (l g |U) 2 )- 1 /2( &/i _ (l g \l h )b g ) and likewise for the 
corresponding creation operator we then have 

H = yigdhsjl - <l 9 |l ft ) 2 {b\b s - b g b\ g ) S z . (15) 

If the electron is in an eigenstate of the ^-operator, the 
effective hamiltonian is exactly the beam-splitter inter- 
action between the two oscillator modes, which for the 
appropriate interaction time will swap not only qubit 
states but, in fact, any quantum states between the g- 
and g-oscillators. With the electron in the spin down 
eigenstate, we have for example 

e~ tHT \i l g g ) = cos0|| l fl a ) - sin0|4 g l g ), (16) 

where 9 = j£l g Qh,y/l — (l g \lh) 2 T. If we want to swap 
an excitation from one oscillator to the other we just set 
= it/2. We denote the corresponding time evolution 
operator Uk . 

From this it is clear that the for a given the time 
needed for the transformation grows as T oc (1 — 
(l g \lh) 2 )~ 2 . To get a sense of how big the overlap might 
be, we can consider the simple case where the potential in 
the dot is 2D harmonic with the ground state ip g (x, y) — 
V'oo^j ?/) & n d excited state iph(x,y) — i()io(x,y), where 
the subscripts counts the number of excitations in the 
x- and y-direction. For this system we get an overlap 
= 3-2 ~ 0.58. 

It is now possible to adress both storage-modes of our 
quantum memory: We first read an arbitrary qubit state 
into the <?-oscillator and we then swap it to the orthogo- 
nal g-oscillator with the pulse sequence described above. 
With the qubit safely stored we can read another qubit 
into the g-oscillator without disturbing the first one, and 
with the swapping mechanism at our disposal we can gain 
random access to any of the two qubits, and using the 
electron-nuclear spin interaction, we can also implement 
quantum gates on the two-bit register. 

V. FIDELITY OF QUANTUM STATE 
TRANSFER 

In this section we investigate how well our pulse se- 
quence Un achieves the ideal time evolution, U^, taking 
into account the higher order terms, that were neglected 
in Eq. (11). The Campbell-Baker-Hausdorff formula is 



exact in the limit of infinitesimal r, and the purpose of 
the current analysis is to assess how fast, and with how 
few steps, we may carry out the total operation without 
significant loss of fidelity. 

We must hence calculate the fidelity between the 
desired and the actually achieved final state, / = 
\(i(!\UiUN\'4')\ 2 , and subsequently average this quantity 
over the relevant initial states of our protocol. Our 
Hamiltonian interaction operators conserve the total 
number of excited spins, and we restrict ourselves here 
to the Hilbert subspace with between zero and two spin 
excitations spanned by 

Uo s O s ),|| 1 9 S )J|0 9 1 § ),|| l g l g )\ || 2 g g ),\i g 2 g ), 

" v ' 

T 

It0 9 9 ),|tl 3 5 ),|t0 s l 9 ). (17) 

We will obtain the exact unitary time evolution in this 
subspace numerically, and we will compare it with the 
desired evolution, applied on states from the the relevant 
input space, i.e., the space of state carrying the four pos- 
sible two-bit states of the spin oscillators. The bracket 
with the symbol S in (17) encompasses this information 
carrying subspace, and we will average the fidelity, / over 
a uniform distribution of input states from S. Since the 
ideal operation restricts the dynamics to the subspace S, 
leakage of excitation to states outside S populated dur- 
ing the exact evolution transformation will automatically 
reduce the fidelity. 

Defining the projection operator, P$ on S , the average 
fidelity can be expressed as the following integral over 
states in S, 

F=(f)= f MP s uiP s U N P s \iP)\ 2 dV, (18) 

where the projection operators are not strictly needed as 
the states are taken from S, and Ul keeps states within 

2 

the subspace. 
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FIG. 4. The infidelity 1 — F is plotted (on a log scale) as a 
function of the number of applied pulses, N, for five different 
values of the overlap (1 9 |1^). 

Using the results in [21] we obtain the explicit expres- 
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F = 



n S {n s + 1) 
/t 



{Tr(AfM t ) + |Tr(M)| 2 }, (19) 



where M = PsUIPsUnPs is a product of projection 

2 

and unitary time evolution matrices, and rig = 4 is the 
dimension of the subspace 5. 

The infidelity, 1 — F, is plotted as a function of the 
number of repeated applications of the pulse sequence 
(10) in FIG. 4, for different values of the overlap between 
the two modes. As we would expect, the infidelity decays 
faster with N when the overlap is small, which reflects 
the fact that it is then easier to address the two modes. 
In the case of a vanishing overlap, we can address the two 
modes independently in a direct manner, but as the over- 
lap increases, it takes longer time and more operations, 
cf., the expression for after Eq. (16). 

Part of the infidelity is due to population escaping the 
subspace S, and ending up with state vector components 
along the other states listed in (17). Since components 
with the electron spin up after the protocol are defi- 
nitely erroneous, we suggest to perform a measurement 
of the electron spin after the swapping operation. With a 
probability that does not exceed the infidelity, plotted in 
FIG. 4, we may find the electron in the spin up state and 
discard the system, while if we find the electron in the 
spin down state, we have effectively removed the small 
electron spin up component from the state, and obtained 
a new normalized state vector belonging to the subspace 
T indicated by the bracket with the same symbol in (17). 

Thus, with a minor reduction of the success probabil- 
ity, we effectively enhance the desired state vector com- 
ponent in S and thus improve the fidelity. Since the 
Campbell-Baker-Hausdorff formula neglects higher order 
terms, which contain the electron spin flip operators, we 
expect that the lowest order error source is associated 
with leakage outside of T, and our heralding may thus 
significantly improve the fidelity of the protocol. 

To determine the average fidelity of the heralded state 
transfer process, we must calculate the integral over all 
initial states from S, of the squared overlap between the 
final (normalized) state projected into T and the target 
state, weighted by the probability of actually finding the 
final state in T [22]. The integral should finally be nor- 
malized by the average probability for the final state to 
be in T: 



Is 



\P T U N P s \ip)fdV 



S s \\P r U N P s m 2 dV 



(20) 



The numerator can be simplified, and using that PsPt 
Ps, we can express the mean conditioned fidelity as 



F r = 



1 Tr(MAft) + |Tr(M)| s 



n s + l Tt(P s U ] N P T U N P s ) ' 
where the operator M is the same as in Eq. (19). 
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FIG. 5. The conditional infidelity 1 — F c is plotted as a func- 
tion of the number of applied pulses, N, and for five different 
values of the overlap (l g |lh). 



We have calculated and plotted the conditional infi- 
delity as a function of the number N of pulse seqeuences 
in FIG. 5 for different values of the overlap (l g |l/j). The 
success probability of the heralding is higher than the 
unconditional fidelity, represented in FIG. 4, and, e.g., 
when F — 99 %, we observe a factor one hundred further 
reduction in the infidelity, associated with the condition- 
ing. 

We have plotted the fidelity as a function of N, but it is 
also of interest to calculate the number of pulse squences 
and the total time, T, required for the protocol at a given 
fidelity level. It is a straightforward exercise to show that 



T 



8ttN 



si g n hy /i - (l g \l h )< 



(22) 



To reach fidelity of F c — 0.999 for two values of the over- 
lap, 0.1 and 0.9, we thus need 



#0.999(0.1) = 19 
#0.999(0.9) = 80 



r~ 2i.9(fi 9 n h )~5 

T~ 67.9(0,^4)- 5. (23) 



As can also be read from the figure, relatively fast swap- 
ping between the oscillator modes is possible with rela- 
tively few operations. In both cases the heralding prob- 
ability of finding the electron in the spin down state is 
close to 95 %. 



VI. DISCUSSION 

In summary, we have shown that ensembles of spins 
can be used as a two- mode storage device. If two orthog- 
onal spin excitation modes can be selectively addressed, 
the operation of the memory is easy, while for a cou- 
pling that addresses non-orthogonal modes, more elabo- 
rate control is needed to identify and subsequently ma- 
nipulate superposition states. 
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We showed that for a wide range of values for the mode 
overlap, which may be implemented with different spatial 
electron wave functions in a quantum dot, it is possible to 
address, swap, and manipulate two qubits of information 
in an ensemble. With more than two choices of spatial 
coupling amplitudes, we imagine that the scheme may be 
readily generalized to more qubits. We further note, that 
our key theoretical component is a beam splitter opera- 
tion, which swaps not only qubit states but also general 
oscillator states between the collective modes, and thus 
provides the possibility to implement controlled opera- 
tion on two-qudit (multi-level) [23-25] and on continu- 
ous variable [2G] states. Our fidelity analysis was carried 
out for qubit degrees of freedom, but we expect that it is 



representative for the unconditioned and the conditioned 
fidelities attainable in qudit and continuous variable sys- 
tems as well. 

Our physical example dealt with the case of an elec- 
tron quantum dot, where the electron spin couples to the 
nuclear spins in the host material within the range of the 
spatial electronic wavefunction. The general problem of 
control of non-ortogonal collective modes may occur in 
different microscopic and mesoscopic systems including 
nuclei in the vicinity of NV centers in diamond, small 
atomic ensembles, and superconducting elements coupled 
to transmission waveguides or to nanomechanical devices, 
and we imagine that multi-mode storage along the lines 
presented here may be pursued in such systems. 
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